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0\ Introduction 

^ ' As it is well known the refraction of the electromagnetic field associated with an electron 

pL| ■ passing through a matter with a uniform velocity originates the Vavilov-Cherenkov and 
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transition radiation. The Cherenkov emission of photons by a charged particle occurs 
whenever the index of refraction n{uj) > 1 (cj is the photon frequency). For the X-ray 



^ I frequencies being higher than atom's ones the refractive index has the general form 

I : "(-) - 1 - 0- m 

^ ' where ujl is the Lengmuir frequency. 
X , ^ 

QJ , But as it was shown for the first time in |I| in the X-ray region the crystal structure 



essentially changes the refraction index for virtual photons emitted under diffraction. This 
leads to the change of generation conditions for bremsstrahlung and Cherenkov radiation. 
As a result, when a charged particle is moving at a constant velocity it can emit quasi- 
Cherenkov X-ray radiation [|1], 0. Moreover, the particle can emit photons not only at 
the angle •& ~ 7~^ (7 is Lorentz factor of the particle) but at the much larger angle as 
well and even in the backward direction. This proved to be the important circumstance 
that allowed to observe Parametric X-ray Radiation (PXR) for the first time |P, ^. The 
observed radiation peaks were shown ||^, |^, ^ to be really a new type of radiation and not 
diffraction of bremsstrahlung radiation. 

In accordance with [2-4] the diffraction of virtual photons in crystal can be described by 
a set of refraction indices n^(u;, k), some of which may appear to be greater than unit. Here 
n{u, k) is the refraction index of crystal for X-ray with photon wave vector k . Particularly, 
in the case of two-wave diffraction the refraction index ni{uj, k) > 1 and n2{uj, k) < 1, and, 
respectively, two waves propagate in crystal - the fast one {n2 < 1) and the slow one 
{ni > 1). For the slow wave the Cherenkov condition can be fulfilled 

1- f3n{uj,k)cos'd = 0, (2) 



where ■& is the angle between photon wave vector k and the particle velocity v (polar angle 
of radiation), (3 = v/c, c is the velocity of light. As a result, this wave can be emitted all 
along the particle pass in crystal. Condition (2) is not valid for the fast wave, and this wave 
can be emitted only at the vacuum- crystal boundary. When the diffraction conditions are 
not fullfilled the well-known ordinary transition X-ray radiation (TXR) is formed. 

PXR is highly directed, quasimonochromatic, tunable, polarized radiation with high 
spectral intensity. These properties as well as the possibility of radiation at the angles 
being much greater than typical radiation angle of a relativistic particle 7"^ make PXR 
potentially useful for wide variety of applications such as materials analysis, microscopy, 
nanofabrication and digital subtraction imaging. Predictions for the spectral and angular 
features of PXR have been corroborated by p, p, [1^, |T^. 

In the above referenced papers the detailed analysis for a diffraction maximum at the 
angle i} » 7^^ has been performed. So in some contributions it has been concluded that 
PXR is emitted at a large angle only |]12|, |n| and it can be described in the kinematic 



approximation by Ter-Mikaelian formula [15] for propagation of radiation in a periodic 



medium. According to [|T4| PXR is the well-known resonance radiation |T5[ . According to 
Q there is a new "PXR of type B" for low energy electrons, which is not related to the 
Cherenkov effect. 

However, it should be noted that according to Baryshevsky and Feranchuk |^, |T^ , each 
photon with the frequency u, emitted at a large angle, corresponds to the photon with 
the same frequency u, emitted at a small angle ■(? < 7^^ near forward direction (particle's 
velocity direction). The frequency of this photon does not depend on particle frequency 
in contrast with resonance radiation photons, frequency of which is proportional to 7^. 
Moreover, the PXR at small angles cannot be described in the kinematic approximation 
by Ter-Mikaelian formula in principle. The kinematic theory cannot describe backward 
Bragg diffraction either. 

Up to now PXR at small angle near the relativistic particle velocity direction has not 
been observed experimentally. PXR observation in the forward diffraction maximum would 
uniquely prove out the quasi- Cherenkov mechanism of PXR. 

In [1^ the expressions for angular distribution of PXR in forward diffraction maximum 



neglecting absorption were derived. For Laue geometry the total number of photons in 



forward maximum was obtained. 

The assumption of low absorption does not allow to take into account the possible total 
reflection of quanta in crystal at Bragg geometry. Because even at negligibly small ordinary 
absorption, the rapidly damped inside the crystal wave arises (so called the internal ex- 
tinction phenomenon |T^). In the present paper PXR for both the forward and backward 
diffraction in Bragg geometry was studied. It was shown that contributions of both the 
slow and the fast waves in radiation intensity are comparable and only their joint account 
allows to describe the experimental results. The comparison of the theory and experiment 
was carried out. The calculations for Si and LiH crystals was performed. 

1 General expressions for PXR emission rate. 

Let a particle moving with a uniform velocity be incident on a crystal plate with the 
thickness L, L << L^r (where Lhj. = (^)~i/2 is the coherent length of bremsstrahlung 
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radiation, q = 9^/4 and 9^ is a root-mean-square angle of multiple scattering of charged 
particles per unit length). This requirement allows us to neglect the multiple scattering of 
particles by atoms of crystal. The theoretical method for case of intense multiple scattering 
can be found in ||^. 

The expression for the differential number of photons with the polarization vector e^ 
emitted in the k direction was obtained in, for example, |ll6| : 



(LjdVt An'^hc^ 



vEi ''^{f{t),uj) exp{—iLut)dt 



(3) 
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where eQ is the particle charge, Ej^ is the solution of homogeneous Maxwell's equation. 
In order to determine the number of quanta emitted by a particle passing through a crystal 
plate one should find the explicit expressions for the solutions Ej^ . The field E^ can 
be found from the relation E J''^ = (Ei )* using the solution E^. , that describes the 

— fe ^ fe ' k 

ordinary problem of photon scattering by a crystal. 

The following set of equations for wave amplitudes in the case of two-beam diffraction 
can be obtained: 
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I-Xl Ei-''-CsX-rEp' = 0. (4) 



Here k^- = k + t , r is the reciprocal lattice vector (|r| = 2TT/d, d is the interplanar 
distance), XojXr'Xir ^^^ ^^^ Fourier components of the crystal susceptibility. It is well 
known that a crystal is described by a periodic susceptibility (see, for example [p!8| ): 

X{.r) =^Xr exp(irf). (5) 

r 

Cs = CsCrs , Cs (eVs) are the unit polarization vectors of the incident and diffracted waves, 
respectively. 

The condition for the linear system (4) to be solvable leads to a dispersion equation 
that determines the possible wave vectors fc in a crystal. It is convenient to present these 
wave vectors in the form k^s = k + n^sN, n^s = -^^tisi where /x = 1, 2; A^ is the unit 
vector normal to the entrance crystalline plate surface and directed inward a crystal. 



^ ^s 



\ [-anPi + Xo(A + 1) ± V[-«bA + Xo(/9i - 1)]' + 4/?iX^xir} , (6) 
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2fcr + r 
«B = p , (7) 

aB is the off-Bragg parameter {as = 0, if the exact Bragg condition of diffraction is 

fulfilled), pi = 7o/7i, 7o = ^7^, 7i = n^rN, n^ = |, n^^ = ||±^. 

The general solution of equation (4) inside a crystal is: 

2 
4~^' (^) = E [e"^M exp(i4sf) + e^A,^ exp(iA^^,,f)] . (8) 

Matching these solutions with the solutions of Maxwell's equation for the vacuum area 
we can find the explicit form of Ej^ (r) throughout the space. It is possible to discriminate 
several types of diffraction geometries, namely, Laue and Bragg schemes are the most well- 
known. 

Let us consider the Bragg case. An electromagnetic wave emitted in the forward di- 
rection as well as the diffracted one emitted by a charged particle and leaving the crys- 
tal through the surface of the particle entrance can be observed. Fig.l demonstrates the 
schemes for forward (a) diffraction maximum and diffraction maximum at large angle equal 
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Figure 1: Schemes of diffraction geometries: a) Forward maximum, b) Diffracted maxi- 
mum. 



to 26b (here 9b - Bragg angle, siiiOb = -^) with respect to particle's velocity direction 
(b) in Bragg geometry. 

Matching the solutions of Maxwell's equations on the crystal surface by use (4), (6) 
and (8), one can obtain the following expression: 
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(9) 



(10) 



(2e2(i)s - Xo) - (2ei(2)s - Xo) exp(-i^(e2(i). - £1(2).)^) ' 

According to (4) in case of two-beam diffraction there are two different refraction indices 
for photons propagating inside a crystal and one of them being able to become greater 
than unit. That is, the superposition of slow and fast waves propagates in a crystal. The 
important thing is that the refraction indices for both waves strongly depend on frequencies 
and propagation directions of photons. 



Substituting (9) in (3) one can obtain the differential number of PXR quanta with 
polarization vector e^ diffracted forwards for the Bragg case: 
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For PXR quanta with the polarization vector eVs diffracted at a large angle relative to 
the charged particle velocity the above formula can be rewritten as follows: 
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2 Transition radiation and PXR 

Let us compare the obtained formulae with the well-known ones for x-ray radiation tran- 
sition in amorphous medium. The TXR formula is known as follows: 
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where ka = k + ^^N is the photon wave vector in amorphous medium. In amorphous 
medium or in a crystal away from the diffraction conditions the only single wave propagates. 
But as a result of diffraction a coherent superposition of several waves propagates in a 
crystal. One can see that (13) is very similar to the expression describing PXR. When 
the frequencies and angles in (11) do not satisfy the diffraction conditions, the formula for 
forward PXR turns to the TXR one (13). At the same time the spectral-angular intensity 
(12) appears to be equal to zero, that is quite obvious since no diffraction at a large angle 
arises. 



3 PXR at small angles with respect to particle's ve- 
locity in Bragg geometry. 

The expression for PXR angular distribution in forward direction can be derived from (11) 
using some assumptions. First of all it should be noted that Xo < and from (6) one 



can obtain that only one root of e^s gives the refractive index n > 1. As a result, the 
difference {u — k^sv) for this /i can be equal to zero and the term in (11) containing this 
difference will grow proportionally to L. On the other hand, the term in (11) containing 
the same difference {u — k^sv) for the other root will not be proportional to L, because this 
difference can never turn to zero. At the first glance it means that the term containing the 
first difference as well as for the Laue case gives the principle contribution to the radiation 
intensity when the particle waylength in crystal Lq grows. But the distinctive feature of 
Bragg case is the possibility of total refiection of quanta that is stipulated by existence of 
the nonhomogeneous wave in a crystal. So, as a result, the both roots of e^s should be 
considered. Therefore, to obtain the PXR angular distribution it is necessary to perform 
the numerical integration of (11) over frequencies in the vicinity of ub- 

Assuming that Lq » I {I = A7^ is the vacuum coherent length, A is wave length 
of the photon) and neglecting the imaginary parts of crystal susceptibilities and possible 
total refiection we can obtain the approximate expression for PXR angular distribution in 
forward direction: 
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We have investigated two geometries for PXR observation in Bragg case in forward 
direction. 

1. As one can see from (14), maximum in angular distribution expression is reached at 
condition 
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(15) 



2c(7-2 - Xo + ^'? 
and the beats over polar angle d may be observed. Obviously it is possible to choose a 

geometry where these effects would be the most essential. Fig.2 demonstrates the experi- 
mental scheme realizing the above condition. 

We have carried out calculations for Si and LiH crystals. Experimental parameters are 
given in table 1. All calculations were performed for azimuth angle if = 7i/2 (this corre- 
sponds to o"-polarization of PXR) and diffraction plane (111). 
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Figure 2: Geometry for observation of beating in PXR angular distribution; 70 
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Figure 3: Spectral- angular distribution of PXR for polar angle -(9 = 2, 09 mrad (azimuth 
angle if = tt/2). 1 - without multiple scattering (MS) taken into account, 2 - with MS. 

Fig. 3 demonstrates spectral-angular distribution of parametric x-radiation for experi- 
mental parameters shown in the 8-th row of table 1. In this geometry PXR is generated 
within frequency range where the coefficient characterizing the efficiency of X-ray reflection 
by the atomic planes becomes nearly equal to unit. In flg.3 the right part of spectrum is 
a flat curve - this is the range of total Bragg reflection. 

The angular distribution for geometry described by the 3-d row of table 2 is shown 
in flg.4. The curve 1 is derived by integration of spectral-angular distribution over the 
frequency range /\uj/ujb = 10^'^ within tuo, corresponding to the center of PXR diffraction 
maximum for each pair {<^,'d). The curve 2 is given by (14), it shows more pronounced 
beating over polar angle than curve 1. The maxima of angular distribution for analytical 



Table 1: Experimental parameters. 



N 


Crystal 


Diffraction plane 


WB,keV 


70 


/3i 


Lq, cm 


Lhr, cm 


Labs, cm 


1 


Si 




30 


0,1211 


-11,57 


8,26-10-2 


9,33-10-3 


0,362 


2 


Si 




30 


0,1177 


-8,43 


8,50-10-2 


9,33-10-3 


0,362 


3 


Si 




30 


0,1142 


-6,55 


8,76-10-2 


9,33-10-3 


0,362 


4 


Si 




30 


0,1055 


-4,03 


9,48-10-2 


9,33-10-3 


0,362 


5 


Si 




30 


0,0882 


-2,02 


1,13-10-1 


9,33-10-3 


0,362 


6 


Si 




30 


0,0659 


-1,00 


1,52-10-1 


9,33-10-3 


0,362 


7 


Si 




30 


0,0349 


-0,36 


2,87-10-1 


9,33-10-3 


0,362 


8 


Si 




20 


0,1636 


-4,85 


3,06-10-3 


1,14-10-2 


0,106 


9 


Si 




20 


0,1636 


-4,85 


3,06-10-2 


1,14-10-2 


0,106 


10 


LiH 




15 


0,3292 


-20,10 


3,04-10-2 


4,60-10-2 


17,4 


11 


LiH 




20 


0,2507 


-25,35 


3,99-10-2 


3,98-10-2 


41,1 


12 


LiH 




30 


0,1707 


-46,71 


5,86-10-2 


3,25-10-2 


139 



formula exceed the maximum values of distribution obtained by integration. This can 
be explained by more exact account of absorption effect while performing integration of 
spectral-angular distribution. Analytical formula was obtained with some assumptions 
(neglecting of imaginary parts of crystal susceptibilities) that leads to some overestimation 
of PXR angular intensity. 

Fig. 5 demonstrates angular distributions of PXR for different values of geometrical 
factor j3i - first seven rows of table 1 (this corresponds to different angles between (111) 
plane and surface of crystal plate). 

2. Another variant of experimental geometry is given by scheme shown in fig. 6. The 
Bragg angle 6b is equal to 90 degrees, this is the geometry realizing normal incidence of 
particles on diffraction plane, the surface of crystal plate is parallel to diffraction planes. 

The calculations for the experimental parameters presented in table 2 were carried out 
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Figure 4: Angular distribution of PXR derived on the base of: 1- integration of spectral- 
angular distribution over frequencies, 2 - analytical formula for angular distribution. 
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Figure 5: PXi? angular distributions for different values of (3i factor (from the left (3i 
-2, 02; -4, 03; -6, 55; 8, 43; -11, 57. 
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Figure 6: Experimental geometry realizing perpendicular incidence. 

The spectral-angular distribution for geometry corresponding to the 3-d row of table 2 
is shown in fig. 7. 

Fig. 8 presents PXR angular distribution corresponding to the 3-d row of table 2. Such 
essential difference between results obtained on the base of analytical formula and that of 
integration of spectral-angular distribution can be explained by the fact that the formula 
(14) does not take into account the nonhomogeneous waves arising at the total reflection 
in Bragg case (in contrast with Laue geometry where the total reflection does not realized) 
(see fig. 7 - the flat part of spectral-angular distribution curve). 

4 PXR at large angles with respect to particle's ve- 
locity. 

As well as for forward PXR, in order to obtain angular distribution for diffracted PXR (see 
fig.l, (b)) taking into account the effect of absorption, possible total reflection of quanta in 
a crystal and contribution of both slow and fast waves in radiation intensity it is necessary 
to perform the numerical integration of spectral-angular distribution (12) over frequencies 
in the vicinity of ujb- We have carried out calculations for the geometry presented in fig. 9. 
Figs. 10 and 11 show spectral-angular and angular distributions for backward diffrac- 
tion. How one can see from Fig. 12, expression (12), taking into account effects of a dynamic 
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Figure 7: Spectral- angular distributions of PXR for polar angle i^ = 12^23mrad (azimuth 
angle) ip = 7r/2 i - without account of MS; 2 - with account of MS. 
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Figure 8: Angular distributions of PXR obtained on the base of: 1 - analytical for- 
mula for angular intensity, 2 - integrating spectral- angular distribution over frequencies 
(Alu/ujb = lO^^y) without taking into account multiple scattering, 3 - integrating spectral- 
angular distribution over frequencies (Alj/ujb = 10^^^ taking into account multiple scat- 
tering. 
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Figure 9: Geometry of backward diffraction. 
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Figure 10: Spectral- Angular distribution for backward diffraction in Si crystal, diffraction 
plane (111), particle's energy 86,9 MeV, "i? = 18,6 mrad (azimuth angle (p = tt/2 ), ujb = 
1, Q'&keV , L = 1, 25 ■ 10~^ cm. 1 - PXR spectral- angular distribution (SAD) for the first 
dispersion branch ji = 1 without the terms proportional to — ^p; 2 - SAD for ji = 2 without 
the terms proportional to — ^; 3 - SAD for the both dispersion branches without the terms 
— ^; 4 ^ SAD for the both dispersion branches without the terms proportional to 
5 - SAD described by the formula (12). 
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Figure 11: Angular distribution for backward diffraction, derived by integration (12) over 
frequencies in the range /\uj/ujb = 10^'^ in silicon crystal for azimuth angle if = 7r/2, 
diffraction plane (111), particle's energy 86,9 MeV, ujb = 1,98 keV, L = 1,25 ■ 10^'^ 
cm. 1- PXR angular distribution for the both dispersion branches without integration the 
terms proportional to — ^; 2- PXR angular distribution with taking into account during 
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integration all terms in expression (12). 
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Figure 12: Angular distribution for backward diffraction, in silicon crystal for azimuth 
angle (p = 7r/2, diffraction plane (400), particle's energy 80,5 MeV, ojb = 4,57 keV, 

results of /|7^/; b - derived by integration (12) over frequencies in 
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Table 2: Experimental parameters. 



N 


Crystal 


Diffraction plane 


WB,keV 


70 


/31 


Lq, cm 


Lbr, cm 


Labs, cm 


1 


Si 


(555) 


9,887 


1,00 


-1,00 


0,05 


1,63-10-2 


1,28-10-2 


2 


Si 


(555) 


9,887 


1,00 


-1,00 


0,01 


1,63-10-2 


1,28-10-2 


3 


Si 


(555) 


9,887 


1,00 


-1,00 


0,005 


1,63-10-2 


1,28-10-2 


4 


LiH 


(333) 


7,876 


1,00 


-1,00 


0,05 


6,35-10-2 


2,52 


5 


LiH 


(333) 


7,876 


1,00 


-1,00 


0,01 


6,35-10-2 


2,52 


6 


LiH 


(333) 


7,876 


1,00 


-1,00 


0,005 


6,35-10-2 


2,52 



diffraction, allows to describe the experiment more precisely than the formula of the kine- 
matic theory [|l5l. 



5 Effect of multiple scattering on PXR. 

In addition to absorption of photons in crystalline plate multiple scattering of charged 
particles by atoms of crystal is the other essential factor also leading to limitation of 
longitudinal dimensions of PXR generation range. So, the PXR characteristics strongly 
depend on relation between crystal thickness Lq and coherent length of bremsstrahlung 
radiation Lbr. When Lq << Lbr the influence of multiple scattering reduces actually to 
a small addition to PXR intensity due to bremsstrahlung radiation mechanism. In the 
opposite case, when Lq > L^r, multiple scattering significantly changes PXR parameters. 
The calculations taking into account multiple scattering effect were carried out on the base 
of expressions given in . 

The relation between L^r and crystal thickness Lq determines the influence of multiple 
scattering on PXR characteristics (the values L^r are contained in the 8-th columns of 
tables 1 and 2) Figs. 3, 7, 8 represents PXR (spectral-angular and angular) distributions 
taking into account multiple scattering influence. All these distributions were obtained for 
thin crystalline plates {L < Lbr)- 
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Figure 13: PXR angular distribution in silicon crystal: 1 - PXR without account of MS; 2 
- PXR with account of MS; 3 - Bremsstrahlung caused by MS. 

In fig. 13 angular distributions of PXR and bremsstrahlung radiation caused by mul- 
tiple scattering for the parameters presented in the 3-d row of table 1 are shown. These 
distributions were obtained by integration of spectral-angular distributions while taking 
multiple scattering into account and without it. This is the case of intense MS, Lq ~ 9 
Lbr, MS leads to strong reduction of PXR intensity. 

Fig. 14 demonstrates dependence of intensity of angular distribution on the thickness of 
crystalline plate. The parameters describing these situations are given in rows 1-3 of table 
2. For plate of thickness Lq = 0, 005 cm MS is not very essential (see fig. 7), for thickness 
of 0,01 cm (Lq ~ Liyj) MS is already noticeable, and in plate with the thickness of 0,05 cm 
MS plays the important role in PXR intensity reduction. 

6 PXR in LiH crystall. 

Use in experiments of LiH crystal enables to realize weak MS even for sufficiently thick 
crystals (see rows 11 and 12 of table 1 and rows 4, 5, 6 of table 2). In case of Bragg 
diffraction for geometry described by fig. 2 it is possible to observe beating over polar angle 
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Figure 14: PXR angular distribution in silicon crystal with account of MS: 1 - L = 0, 005 
cm, 2 - L = 0,01 cm, 3 - L = 0, 05 cm. 

in angular distributions measurements. The experimental geometries more advantageous 
for PXR observation are marked in table 1 by (*). For these experimental schemes the 
comparison between angular intensities of PXR and TXR was carried out {Auj/ub = 10"'*). 
Figure 15 demonstrates angular distributions of PXR and TXR for LiH crystal. 

For Bragg diffraction in case of normal incidence the geometries preferable for PXR 
observation possibility are asterisked at margins of table 2. 

Fig. 16 demonstrates comparison of PXR and TXR angular intensities for situation 
described by the 4-th row of table 2. 

Conclusion. 

The only precision account of all dynamic effects (total Bragg reflection, interference of 
slow and fast waves or, by the other words, taking into consideration of all terms in (12), 
allows to obtain the real picture of diffraction radiation in case of Bragg geometry. The 
contributions of both the slow and the fast waves in radiation intensity for Bragg geometry 
are comparable and only their joint account allows to describe the experimental results. 
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Figure 15: Angular distribution in LiH crystal for geometry described by fig .2: 1 - PXR, 
2 - TXR. 
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Figure 16: Angular distribution in LiH crystal, case of normal incidence geometry: 1 
PXR, 2 - TXR. COB = 7, 88 keV, Lq = 0, 01 cm. 
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